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Reverse— domain superconductivity in superconductor— ferromagnet hybrids: 
effect of a vortex— free channel on the symmetry of I — V characteristics 
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We demonstrate experimentally that the presence of a single domain wall in an underlying fer- 
romagnetic BaFei2 0ig substrate can induce a considerable asymmetry in the current (7) - voltage 
(V) characteristics of a superconducting Al bridge. The observed diode-like effect, i.e. polarity- 
dependent critical current, is associated with the formation of a vortex-free channel inside the 
superconducting area which increases the total current flowing through the superconducting bridge 
without dissipation. The vortex-free region appears only for a certain sign of the injected current 
and for a limited range of the external magnetic field. 



PACS numbers: 74.25.F- 74.25.Dw 74.78.Fk 74.78.Na 

The development of material deposition techniques 
and lithographic methods have made it possible to fab- 
ricate superconductor-ferromagnet (S/F) hybrid struc- 
tures with controlled arrangements of ferromagnetic 
layers/elementsji - — These flux- and exchange-coupled 
S/F hybrids^— are of fundamental interest for investi- 
gations of nontrivial interactions between superconduc- 
tivity and nonuniform distributions of magnetization. In 
addition S/F hybrids seem to be potential candidates for 
the development of tunable elements of superconducting 
electronics^ 

It is known that a nonuniform magnetic field can mod- 
ify the conditions for the appearance of superconduc- 
tivity due to the effect of a local field compensation^ 
In flux-coupled S/F bilaycrs^ the formation of localized 
superconductivity results in either domain-wall super- 
conductivity (DWS) or reverse-domain superconductiv- 
ity (RDS), when superconductivity occurs, respectively, 
above magnetic domain walls or above magnetic domains 
of opposite polarity with respect to the orientation of an 
external magnetic field H ext (see review^ and references 
therein) . The appearance of localized superconductivity 
(DWS and RDS) becomes possible if the amplitude of the 
nonuniform field, Bq, is comparable or exceeds the up- 
per critical field, H C 2, of the superconducting material, 
which was confirmed experimentally for various planar 
S/F structures £"-12 

A present challenge, associated with these S/F hybrids, 
is the direct investigation of the transport properties of 
superconducting channels that arc induced by stray mag- 
netic fields. Indeed, this problem seems to be crucial for 
any practical applications exploring the effect of localized 
superconductivity and guided vortex motion in tunable 
magnetic landscapes. Parallel magnetic domains in thick 
permalloy films were found to lead to a preferential vor- 
tex motion and a giant anisotropy of the critical currents 
in S films and crystals,— ~— even though the amplitude 
of the nonuniform field appears to be insufficient for the 
formation of localized superconductivity (Bq/H C 2 < 1 at 




(b) 




FIG. 1: (color online) (a) Combined magnetic and atomic 
force microscopy images of the hybrid sample. The dark and 
bright stripes, corresponding to the different ferromagnetic 
domains, are shown only within the cross-shaped supercon- 
ducting bridge. White dashed lines depict the position of the 
domain walls. All elements of the electrical circuit is shown 
schematically, (b) The model S/F structure. 



low temperatures). 

It is important to note that the magnetic field induced 
by parallel magnetic domains in BaFci20i9 is rather high 
and thus suitable for RDS in superconducting Al films 
since Bq/H C 2 > 2.5 for all temperatures^ Using such 
S/F bilayers with well defined localized superconducting 
channels in a normal-metal matrix, we continue our pre- 
vious studyH with the aim to test the potential of the 
superconducting channels to carry current. In this Letter 
we focus on the measurements of the current (I) - volt- 
age (V) characteristics of a S/F bilayer along a single 
domain wall in a ferromagnetic substrate as a function of 
H ext and analyze the dependence of the critical current 
I c on H ext for different signs of the bias current. 

Our sample consists of a ferromagnetic crystal 
BaFei20ig with a thin-film superconducting Al bridge 
grown on top. Since the ferromagnetic and supercon- 
ducting parts were electrically isolated by a 5 nm Si 
buffer layer, the interaction between these parts was 
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FIG. 2: (color online) The z— component of the field in- 
duced by the ferromagnetic domains, measured by a scanning 
Hall probe microscope along a line perpendicular to the do- 
main walls at T = 72 K, H ex t = and at height 400 nm. 
The dashed line corresponds to 77/ = 77o arctan (y/L) 
(77,, = 115 Oe, L = 1.5 /im). 



purely magnetostatic. Being cut along the proper crys- 
tallographic direction, the polished crystal BaFe^Oig 
exhibits a stripe-type domain structure with dominant 
in-plane magnetization^ The location of the domain 
walls was determined by magnetic force microscopy, prior 
to the preparation of the superconducting bridge. The 
cross-shaped Al microbridge (30 /xm wide and 50 nm 
thick) was fabricated by e-beam lithography, molecular 
beam epitaxy and lift-off etching [Fig. QJa)]. A similar 
structure was used in RefJ^ for the observation of the 
anisotropy of the electrical resistance in this S/F bilayer. 

The profile of the perpendicular z— component of the 
nonuniform magnetic field is shown in Fig. [2] Its am- 
plitude, even though measured at a rather large distance 
(400 nm) from the surface, is close to the H c i value for 
Al films (H c2 ~ 2 • 10 2 Oe at T = 0). Since the field 
amplitude inside the superconducting film, Bq, will ex- 
ceed H C 2 at all temperatures, superconductivity appears 
in this S/F system only at H ext ^ as reverse-domain 
superconductivity above positively (negatively) magne- 
tized domains at H ext < (H ext > 0), respectively. 

Figure [3] shows the typical I — V characteristics mea- 
sured at T = 0.5 K (T/T c0 ~ 0.34). Depending on the 
Hext value, there are three different cases: 

(i) symmetric normal-type I — V dependence with al- 
most constant slope dV/dl (not shown here); 

(ii) symmetric I — V dependence with non-zero critical 
current (curves labelled 450 Oe and 470 Oe). This is 
realized in a rather wide H ext range corresponding to 
the reverse-domain superconductivity; 

iii) asymmetric hysteretic I — V dependence (curves 

Here 
for the 

ascending branches of the I — V curves both for positive 
(+) and negative (— ) polarities of the transport current. 
This type oi I — V characteristics was found only in the 
close vicinity of the compensation field (|i? ex t| — Bq). 



labelled 490 Oe and 510 Oe) with l| +) + A } 
we introduce the critical currents D and D 
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FIG. 3: (color online) Typical I — V dependencies mea- 
sured along the domain wall for different H ex t values. Both 
branches for 7 > and 7 < were measured starting from 
7 = 0. 
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FIG. 4: (color online) (a) Critical currents 7c as a function 
of H ex t at T — 0.5 K. The shaded areas correspond to the H ex t 
range where the decrease of the resistance from its normal 
value (i.e. the RDS regime) was detected. Arrows indicate 
the compensation field |77 ea; t| = Bo- (b) Comparison between 
experiment (circles) and theory (solid lines). 



are summarized in Fig. @Ia). The relationship between 
Ic and Ic depends both on the absolute value of H ext 
and its sign: 4' > 4 at H ext < and vice versa. 
Thus, the most important finding of this paper is the 
field-induced change of the symmetry of the I — V char- 
acteristics. The fact that the transmission capacity of 
the superconducting channel formed in the non-uniform 
magnetic field, can be strongly dependent on the polar- 
ity of the transport current may look rather unusual and 
counter- intuitive. 

To clarify the physical origin of the difference between 



(a) H ext ~B o J I <0 



(d) H cxt ~B J x >Q 
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FIG. 5: (color online) Static vortex and current patterns 
in the S bridge in the presence of the non-uniform field 
Hf = Ho arctan (y/L) , corresponding to the maximum of 
the flowing current, at different H ex t values. V and AV 
stand for vortices and antivortices, jdep is the depairing cur- 
rent density' . The right part of the superconducting bridge 
(0 < y < w/2) is in the normal state at H ex t > and j = 
and n — there. 



7c and Ic , we analyze the transport properties of 
the generic S/F bilaycr within the London model. Wc 
consider a thin superconducting bridge (width w) in the 
field of a single straight domain wall positioned in the 
center of the bridge (y = 0) and aligned parallel to the 
x— axis [Fig.QJb)]. We assume that all parameters of the 
resulting current /vortex structures depend only on the 
transverse y— coordinate. The current density j(y) and 
vortex density n(y) in the presence of a nonuniform field 
meet the Maxwell-London equatio n 19 ' 20 



to/2 

4^d L + d j j(y')dy' 
c dy c J y 1 - y 

-w/2 



H ext + H f (y)-$ n(y), 



where Hf(y) = Hq arctan(y/L) is the field induced by 
the domain wall, A is the London penetration depth, d 
is the thickness of the superconductor, $ = ^hc/ e is 
the flux quantum. To determine the critical current I c 



of the superconducting strip, corresponding to the tran- 
sition between the state with motionless vortices to the 
flux motion regime, we apply the following conditions: 

(i) the maximum of the current density in the vortex- 
free region should be equal to the critical current density 
j s , which defines the threshold value for the nucleation 
of vortices and antivortices inside the superconductor or 
at its edges; 

(ii) in order to guarantee a flux motion regime at 
I = I c , j(y) should be equal to the depinning current 
density j p in the area where n(y) ^ 0, with j p = 
Jp0 /(l + \B Z \/B p ) 
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is the local magnetic field. As a result, in the vortex-free 
region, n(y) = 0, the current density can be larger than 

Jp- _ 

(iii) the profile J*(y), which satisfies both conditions (i) 
and (ii), allows us to define the critical current as follows 

I c = f™ w 2 /2 3*(y)dy- 

In our calculations we use the parameters typical for 
our system: w = 30 /im, d = 50 nm,A = 150 nm, 
Bo = 520 Oe, B p = 30 Oc, H = 331 Oe, j p0 = 0.14 j dep , 
j s = 0.55 jdep, where jdep is the depairing current den- 
sity of Al at low temperatures^ Our choice for the pa- 
rameter L = 0.35 /im seems to be reasonable since the 
width of the transient area in the i? z -distribution inside 
the superconducting film (at h < 50 nm) can be sub- 
stantially smaller than that measured at large distances 
(L = 1.5 /Urn at h = 400 nm). Since there is rather 
good agreement between the experimental data and the 

calculated dependencies Ic (H ex t) [Fig. 1Kb)], we can 
interpret the asymmetry of the transmission capacity of 
the superconducting channels as follows. 

Provided H ex t — Bo, the magnetic field is effectively 
compensated in the left part of the bridge, while the right 
part will be switched to the normal state. Since the local 
field H z = H ext +Hf changes its sign inside the supercon- 
ducting area (in the case H ex t < Bo), the stable vortex 
structure should generally consist of vortices and antivor- 
tices. However the Lorentz force F^ = c _1 [j x <fr ] act- 
ing on a vortex depends on the direction of the transport 
current j, therefore the resulting vortex pattern, corre- 
sponding to the non-dissipative current flow, may be de- 
pendent on the sign of I x . Indeed, an injection of the neg- 
ative bias current (I x < 0) forces vortices (antivortices) 
to move to the right (left), resulting in a vortex-free chan- 
nel inside the RDS area [Fig.[5](al)-(cl)]. The exact po- 
sition a and the width 5 of such channel depends on H ext . 
Since there are no vortices in a certain area, one can ap- 
ply a larger current through such a vortex-free channel 



without loosing energy and the excess current due to this 
effect can be roughly estimated as S x (j s — j p0 ). How- 
ever, the vortex-free area is absent for a bias current of 
opposite polarity (I x > 0), since vortices and antivor- 
tices move in counter directions and annihilate at the 
point of zero magnetic field [Fig. [5] (a2)-(c2)]. In this 
case the current density cannot be larger than j p o and 

I c < I c . In the under-compensated regime, when 
the absolute H ex t value is substantially less than Bq, the 
gradient of the local field (dH z /dy) a increases rapidly as 
H ext decreases [compare Fig.[S](al)-(cl)], 8 — > 0, and the 
diode effect vanishes. In the over-compensated regime 
(H ext > B ) the vortex-free region positioned near the 
left edge of the bridge becomes very narrow (S — > 0) and, 
as a consequence, the excess current goes to zero and the 
symmetry of the I — V characteristics is restored. Obvi- 
ously that for H ext < we have the same physics, but the 
excess current corresponds to the opposite case (/ > 0) 



and therefore Ic > Ic ■ All these conclusions arc in 
agreement with our experimental observations [Fig. 0] . 

Summarizing, we showed that a nonuniform field can 
cause a pronounced asymmetry of the I — V curves of 
a superconducting bridge provided \H ext \ ~ Bq. The 
difference in the critical currents can be attributed to 
a removal of vortices from the inner part of the super- 
conducting bridge under the action of the Lorentz force. 
Such a vortex-free channel, forming only for a certain 
polarity of the injected current, is able to carry extra 
current without dissipation and thus prevents the super- 
conducting bridge from switching to the normal state. 
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